Perturbative Renormalizability of Chiral Two Pion Exchange 
in Nucleon-Nucleon Scattering 
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We study the perturbative renormalizability of chiral two pion exchange for the singlet and 
triplet channels within effective field theory, provided that the one pion exchange piece of the 
interaction has been fully iterated. We determine the number of counterterms/subtractions needed 
in order to obtain finite results when the cut-off is removed, resulting in three counterterms for the 
singlet channel and six for the triplet. The results show that perturbative chiral two pion exchange 
reproduce the data up to a center-of-mass momentum of k ~ 200 — 300 MeV in the singlet channel 
and k ~ 300 - 400 MeV in the triplet. 
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I. INTRODUCTION 

The effective field theory formulation of nuclear 
forces [l|-m tries to exploit in a systematic manner the 
separation of scales between pion physics, which is known 
to dominate at large distances, and short range physics 
in the two-nucleon system. In Weinberg's original pro- 
posal [a, |6| the chiral nucleon-nucleon potential is orga- 
nized as a power expansion (or counting) in terms of Q 

VnnW = V {0 \r) + vW(r) + V^(r) + 0(Q 4 ) , (1) 

where Q represents the low-energy scales of the sys- 
tem, usually the momentum p of the nucleons and the 
pion mass m w . The potential is then inserted into the 
Schrodinger or Lippmann-Schwinger equation in order to 
obtain theoretical predictions |7H32j . This prescription is 
usually referred to as Weinberg's counting. 

The resulting chiral potentials turn out to be singular, 
behaving at order Q v as l/r 3+lJ in coordinate space for 
short enough distances (m w r C 1). Therefore they need 
to be regularized in order to obtain well-defined results, 
usually by introducing a cut-off in the computations plus 
the necessary number of counterterms which ensure the 
renormalizability of the scattering amplitude Q- This 
has been found to be in contradiction with Weinberg's 
power counting, where the corresponding counterterms, 
determined by naive dimensional analysis, are not able 
to render the theory renormalizable [l9l-[21, 26] (or gener- 
ate chiral inconsistencies [§] , prompting the KSW count- 
ing [33|, Hi| ) . Consequently one is forced to make a dc- 
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cision: either to follow an a priori power counting or 
require renormalizability. 

The direct and practical choice is to follow Weinberg's 
original counting unaltered, leading to a framework ami- 
cable with large numerical computations, which demysti- 
fies nuclear forces and enjoys an undisputed phenomeno- 
logical success [HI, . The price to pay is that the cut- 
off must be fine tuned, lying inside a narrow window, a 
situation which we regard as unsatisfactory from a theo- 
retical viewpoint. Recentl y, b ased on the renormalization 
philosophy of Lepage [H,l3B|, there has been interesting 
attempts to justify this particular approach [37L IHH • 

On the contrary, if one strives for a more robust the- 
oretical foundation, one should be able to achieve cut- 
off independence. The results from non-perturbative 
renormalization in the case of singular interactions 
I2H Ull I39H41I can be summarized as follows: one coun- 
terterm is needed to renormalize a channel where the po- 
tential is attractive and singular, while channels where 
the potential is singular and repulsive become insensi- 
tive to counterterms. The first condition can lead to 
an alarming loss of predictive power, as already at lead- 
ing order (LO) there is an infinite number of attractive 
singular channels. The solution proposed in Ref. [l9| is 
to treat all partial waves with sufficiently high angular 
momentum perturbatively at LO, a procedure which is 
supported by the analysis of Ref. J42[ q The second con- 
dition is particularly problematic: in the triplet channel 
the potential is attractive at LO but becomes repulsive 
at next-to-leading order (NLO), resulting in an unbound 
deuteron at this order when the cut-off is removed pol |. 
As there is no way to predict what the sign of the inter- 
action will be at higher orders, this represents a contin- 
uous threat to the non-perturbative renormalizability of 
the chiral potentials. In addition, there exists the risk 
that non-perturbative renormalization of the subleading 
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pieces of the potential may lead to incompatibilities with 
the chiral expansion [38[. The previous issues can be 
avoided with the perturbative treatment of the higher 
order pieces of the potential, which respects power count- 
ing and renormalizability independently on whether the 
subleading contributions are repulsive or attractive. The 
problem is how to construct such a perturbation theory. 

The purpose of this paper is to investigate the con- 
ditions under which perturbative chiral two pion ex- 
change (TPE) can be renormalizcd in order to extend 
the power counting proposal of Nogga, Timmermans and 
van Kolck 19[ to subleading orders. In the spirit of 
Refs. [H|,|43|, we use renormalizability as a guide to iden- 
tify the required short distance operators. The technical 
meaning of renormalizability depends on whether we are 
in a perturbative or non-perturbative context. By per- 
turbative renormalizability we refer to the elimination 
of all negative (positive) powers of the coordinate (mo- 
mentum) space cut-off in the observables. In contrast, 
non-perturbative renormalizability deals with ambigui- 
ties instead of divergences: the scattering amplitude of 
an attractive singular interaction is finite but non-unique 
and requires the inclusion of a counterterm for fixing the 
solution [2(1 Hi], H^] . A particularly straightforward man- 
ner to fulfill the renormalization program is to study the 
cut-off dependence of the amplitudes when the cut-off is 
removed, as exemplified in Ref. [19]. This should not 
be interpreted however as the necessity of eliminating 
the cut-off in the computations: after the renormaliza- 
tion process, the residual cut-off dependence of the am- 
plitudes is in principle a higher order effect, provided the 
cut-off lies within a sensible range. 

The perturbative techniques in this paper are directly 
based on those sketched in Ref. [2(J, and are equiva- 
lent to the momentum space perturbative methods de- 
veloped in Ref. [43(. Here we use renormalized distorted 
wave Born Approximation (DWBA) with the aim of con- 
structing phase shifts. Complementarily, the approach of 
Refs. [44T447| employs the DWBA techniques for "decon- 
structing" the phenomenological phase shifts, that is, for 
extracting the corresponding short range physics once the 
long range pion effects have been removed and checking 
whether this short range interaction is consistent with the 
specific power counting under consideration, be it either 
Weinberg [44] or Nogga, Timmermans, van Kolck (45l — 
H3] • Of particular interest is the recent deconstruction of 
the 1 S'o singlet channel [46| which advances some of the 
results and conclusions of the present work. The present 
approach differs however from the finite cut-off pertur- 
bative set-up of Ref. [48[ , in which not all the operators 
needed to obtain a renormalized results are included as 
the previous work concentrates on analyzing the Wein- 
berg counting. 

The paper is organized as follows: in Section |TT] we 
study the perturbative renormalizability of the So sin- 
glet channel and determine the cut-off and momentum 
range for which an acceptable description of the data 
is obtained. We extend the previous results to the the 



3 Si — 3 Di triplet channel in Section [TTT1 The role of the 
cut-off within the present approach is analyzed in Section 
ITVl and the relation with other approache s, p articularly 
the renormalization group analysis of Ref. [42] , is consid- 
ered in detail in Section [V] Finally, we briefly summarize 
our results in I VII The technical details of the perturba- 
tive treatment of chiral TPE are explained in Appendices 
IA1 and [Bj Some of the 1 So singlet results from this paper 
have been advanced in Ref. 49]. 



II. SINGLET CHANNEL 

The present perturbative treatment of chiral TPE is 
based on distorted wave Born approximation. For sim- 
plicity, we will only consider in detail the singlet case. 
We can express the phase shifts as the following series 

S(k; r c ) = *C°> (k; r c ) + 6& (fc; r c ) + 5^ (k; r c ) + 0(Q 4 ) , 

(2) 

which is ordered according to the counting of the finite- 
range piece of the potential 0. That is, power counting is 
now manifest in the amplitudes. The LO phase shift 5^ 
is computed non-perturbatively (and includes one coun- 
terterm0), while 5^ and 5^ are computed in first order 
perturbation theory 0. The corresponding expression for 
the perturbative phase shifts is (see Appendix |A"|) 



#">(fc;r c ) 2^ 



sin 2 (5(°) 



= -±A^\k ] r c ) 2 I^(k-,r e ), (3) 



where v = 2, 3 and the perturbative integral l[ V g is de- 
fined as 

/*°° 2 
li V l{k;r c ) = / drV^{r)uf (r) . (4) 



In the previous formulae /i is the reduced mass, k the 
center of mass momentum, is a normalization fac- 
tor, which is taken to be unity at k = 0, and is 
the LO reduced wave function in an energy independent 
normalization at the origin (or at the cut-off radius r c 
if we are using a finite cut-off). The asymptotic nor- 
malization of it^ is determined by A^°\k)u^(k) — > 
sin (kr + sin <j(°) for r — > 00. 

As can be easily checked, the perturbative integral di- 
verges as 1/ r" +2 as a consequence of the short distance 
behaviour of the reduced wave function (r) ~ 1 and 
the potential V {v \r) ~ l/r" +3 . The diver gences can be 



3 Strictly speaking, the leading order piece is of order Q~ 1 and 
not Q°. However, for keeping the notation simpler, we have just 
followed Eq. JTJ). 

4 Note that we are not considering here chiral symmetry breaking 
terms separately. 

5 Second order perturbation theory is not needed as the iteration 
of the NLO potential is of order Q 5 . 
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cured by making the adequate subtractions. Due to the 
energy-independent normalization of u 1 ^ 1 at the origin, 
the terms in the k 2 expansion of uf* = J2 n u 2n^ 2n are 
progressively less singular, with ~ r 2n for r — y 0. Ex- 
panding the previous integrals in terms of k 2 for v = 2, 3 
we have 

(fc; r c ) = /W(r c ) + fc 2 /^(r c ) + fc 4 J<">(r e ) 

+ ^ >H (A;rc), (5) 

where 7q 2,4 are the divergent pieces of the integral and 
R is the regular piece, as can be trivially checked. 
Therefore three subtractions or counterterms are needed 
in order to renormalize the perturbative results in the 
singlet. The specific method employed is not important. 
Here we modify the perturbative integral by adding three 
free parameters which are to be fitted to the data 

r c ) = \t ] + X^k 2 + fc 4 + (fc; r c ) . (6) 

By assuming the short range physics to be parametrized 
by an energy dependent delta-shell potential of the type 

Vc V) r c ) = ^ E C *n (^ k2n (7) 
c n 

we can easily relate the X^J parameters to the coun- 
terterms by 

^ = ^ctt(r c )uP\r c ). (8) 

Equivalently, if one chooses to work in the momentum 
space formulation of Ref. [HJ, one could include the 

contact potential (pl^V) = Cq ] + C^V + p' 2 ) + 
C^\p A +p' 4 ). In either case, the first free parameter, 

Xq (Cg ), is only used to absorb the k = divergence 
of the perturbative integral while its finite piece is redun- 
dant, as it only affects the zero energy behaviour of the 
phase shifts which has already been fixed at LO, meaning 
that we need to fix two additional observables, for exam- 
ple the effective range tq and the shape parameter vi , in 
order to determine the NLO/N 2 LO results. The number 
of counterterms agrees with the corresponding one pre- 
dicted in the RG A of Ref. [42| , where the power counting 
resulting from treating one pion exchange (OPE) non- 
perturbatively was analyzed in detail, and with the re- 
lated deconstruction of Ref. [46[ , in which the short range 
physics for the singlet channel is determined by removing 
the non-perturbative OPE and perturbative TPE effects 
from the phenomenological phase shifts. Note that in 
Ref. |20j an incorrect number of counterterms was deter- 
mined due to an improper normalization. 

The results for the singlet 1 So channel are shown in 
Fig. ©. Following [H E3, we take f« = 92.4 MeV, 
= 138.03 MeV, g A = 1.26, and d ls = -0.97 GeV 2 . 
For the chiral couplings we employ the customary val- 
ues ci = -0.81 GeV~\ c 3 = -3.40 GeV" 1 and c 4 = 
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FIG. 1. (Color online) Phase shifts for the 1 So chan- 
nel with non-perturbative OPE and perturbative TPE. The 
non-perturbative OPE computation contains one countert- 
erm which is determined by fixing the 1 S'o scattering length, 
ao,s = — 23.74 fm, while the perturbative TPE computation 
contains a correction to the LO counterterm plus two ad- 
ditional counterterms which are used to fit the Nijmegen II 
phase shifts [501 ] (equivalent to the Nijmegen PWA [511]) in 
the range k — 40 — 160 MeV. The error bands are gener- 
ated varying the cut-off within the 0.6 — 0.9 fm range. The 
light blue band represents the N 2 LO results from the stan- 
dard Weinberg approach of Ref. [lTj ■ The dashed dark blue 
line represents the N 2 LO results for r c = 0.1 fm. 

3.40 GeV -1 , which are compatible with the determina- 
tion of Ref. [H3 ■ The potential is taken from Ref. [53| . 
As can be seen, the results reproduce the 1 5'o phase shifts 
up to k ~ 200 — 300 MeV, depending on the value of the 
cut-off. If the cut-off is small (r c = 0.1 fm), the perturba- 
tive treatment of the subleading pieces of the interaction 
starts to fail already at k ~ 200 MeV, as a consequence 
of the relative weakness of OPE with respect to the enor- 
mous strength of TPE at short distances. The previous 
problems can be circumvented by using cut-offs of the or- 
der of r c ~ l/2m,r (0.7 fm), which are small enough as to 
guarantee the correct inclusion of the TPE tail. In par- 
ticular we employ r c = 0.6 - 0.9 fm@ a range for which 
perturbative TPE calculations compete well with non- 
perturbative ones in the Weinberg counting at the same 
order [l3|, [TtJ , though perturbative TPE is slightly less 
predictive due to the additional counterterm. On a differ- 
ent ground it should be noticed that OPE is perturbative 
in the singlet [H, [H, HH , even if iterated [55| , suggest- 
ing that the previous results could be reinterpreted as an 
N 3 LO/N 4 LO computation in the KSW counting |33l.l3H. 

The failure of perturbative subleading TPE at r c = 
0.1 fm raises interesting questions regarding the adequacy 
of the present power counting scheme and the role of chi- 
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ous configuration space cut-off range is approximately equivalent 
to a momentum space (sharp) cut-off of A ~ 350 — 500 MeV. 
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ral TPE. Of course, the technical reasons why perturba- 
tion theory fails already at k ~ 200 fm for small cut-offs 
are clear: OPE does not provide enough long range dis- 
tortion as to avoid higher momentum waves to probe the 
van der Waals component of TPE, as has been discussed 
for example in Ref. 46] . This component originates from 
the behaviour of subleading TPE, which in the singlet 
channel can be schematically written as [531 ] 

2/x^pe 3) W = e~ 2m * r E ^m n r) n , (9) 

T 71=0 

where the a n 's are dimensionless parameters with <iq = 1 
and i?6 is a length scale related with the strength of TPE 
at short distances, which varies between Rq = 1.6 — 1.8 fm 
for typical values of the chiral couplings. The previous 
form implies that the chiral van der Waals component 
of subleading TPE should start to become apparent at 
distances below r < l/2m w ~ 0.7 fm. This figure is 
supported by several renormalized non-perturbative TPE 
computations in the singlet [2^, [2l|, [24|, which usually 
reach cut-off independence at distances around or be- 
low 0.5 fm, signalling the onset of chiral van der Waals 
forces. For such cut-off radii the perturbative treatment 
of TPE generates terms like kR^ and m^Re, which, tak- 
ing into account the size of i? 6 , might cause the pertur- 
bative series to eventually diverge. The most consistent 
and straightforward solution to this problem is the use 
large enough cut-offs (r c > 0.5 fm) in order to avoid the 
conjectured breakdown of the perturbative series. The 
alternative solution, which will not be considered in the 
present work, is the iteration of chiral TPE or at least 
some parts of it [46| . Although interesting, this proposal 
seems difficult to harmonize within the EFT framework 
as it requires (i) to justify the promotion of an order Q 3 
interaction to order Q~ l and (ii) the existence of a cut- 
off window for which subleading TPE dominates but the 
higher order corrections are still small compared to this 
contribution. 

The employed cut-off window, r c — 0.6 — 0.9 fm, repre- 
sents a compromise between the requirements of the sin- 
glet and triplet channels. The optimum value of the cut- 
off in the singlet lies in the vicinity of r c = 0.9 — 1.0 fm, a 
range for which the description of the triplet phases starts 
to worsen. This cut-off window may look soft, but it is 
not: the first deeply bound state (i.e. the first zero of the 
k = wave function) for the N 2 LO potential happens at 
r c = 0.70 fm, meaning that the lower range of the present 
cut-off window is already beyond what can be reached in 
the Weinberg scheme. It is interesting to notice that the 
previous cut-off range is similar to the radii at which 
most potential models of the NN interaction (5(J HH, 
have their minima, usually at r ~ 0.8 — 0.9 fm. The 
mimima mark the distance at which the short range re- 
pulsion starts to overcome the long range attraction, and 
consequently can be understood as the separation point 
between short (r < 0.5 fm) and long range (r > 1.0 fm) 
physics. In this sense, the cut-off is to be interpreted as a 



separation scale, as has been proposed within the context 
of RGA [H, [H |B3, rather than as a hard scale [H S|. 

In the calculations of Fig. (JTJ) we also interpret the cut- 
off variation of the results as the error band of the theory. 
The previous is a sensible prospect in the sense that we 
expect the cut-off dependence of the scattering ampli- 
tudes to be a higher order effect. However, if the cut-off 
variation is to be understood as an error band, the size of 
the band should decrease at each new order to reflect the 
convergence properties of the theory. Paradoxically the 
N 2 LO band is bigger than the NLO one, a worrisome sit- 
uation which does not necessarily mean that we should 
abandon the previous interpretation. In fact the same 
happens in the Weinberg counting, as illustrated by the 
singlet channel results of Ref. [ljj • The explanation is to 
be found in the surprisingly large size of the C3 and C4 chi- 
ral couplings, which causes the subleading TPE contribu- 
tion to the chiral potential to be substantially bigger than 
the corresponding one from leading TPE. This is due to 
the large contributions from the A resonance to the chi- 
ral couplings [6(|, c 3 ,a = -2c 4i A = —ih\/9A, with A 
the nucleon-delta mass splitting and Ka the nNA axial 
coupling, ranging from —1.7 to — 2.7GeV~ 1 depending 
on the value of Ha El In this sense, the increased size 
of the N 2 LO error bands is just a reflection of the un- 
expected contribution from this low energy scale. The 
explicit inclusion of the A resonance in the NN chiral po- 
tential, a theme which has been recurrently considered in 
the literature [|, [Q, [H, HUll , is presumed to solve the 
current issue with the error bands (see also the related 
discussion of Ref. [46j]). This prospect does not appear 
to be unreasonable in view of the perturbative peripheral 
wave N 2 LO-A results of Ref. [22[ and the related non- 
perturbative central and peripheral wave calculations of 
Refs. [24], HH, all of which indicate an enhancement in 
the convergence rate of the phase shifts as compared to 
the A-less theory. 



III. TRIPLET CHANNEL 

In the case of the 3 Si — 3 Di channel the perturbative 
analysis is analogous to the previous one for the 1 So chan- 
nel, but more cumbersome due to the presence of coupled 
channels and the singular behaviour of the tensor piece 
of the LO potential in the triplet channels. The details 
of such analysis are shown in Appendix [B] but the essen- 
tial point is that the inverse power law behaviour of the 
OPE tensor force (~ 1/r 3 ) softens the perturbative inte- 
grals and reduce the necessary number of counterterms 
per phase. In fact we have that the s- and d-wave wave 
functions behave as tt£ , wf* ~ r 3 / 4 near the origin [6l| , 



7 Taking h A between 1.08 and the SU(4) value 1.34, see [53 ■ The 
values for the chiral couplings once the A has been included can 
also be consulted in Ref. 12211 . 
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and that each subtraction adds an r 5 / 2 suppression to the 
short distance behaviour 0. This translates into two sub- 
tractions for each of the three phases in the 3 Si — 3 Di 
channel (Sa Sl ,ei, 63 Dl ), meaning that we end up with six 
counterterms at NLO/N 2 LO in agreement with Rcf. 201. 
That is, the scattering amplitude can be completely de- 
termined using six pieces of data, for example the value 
of the three phase shifts at two different momenta. 

The results are shown in Fig. ([2]). Perturbative TPE 
provides a good description of the 3 Si and 3 £) 1 phases 
and the t\ mixing angle up to moderately high momenta, 
around k ~ 300 — 400 MeV, although it should be noted 
that the results are quite sensitive to the choice of chi- 
ral couplings, due to the linear dependence generated 
by treating chiral TPE perturbatively. Contrary to the 
singlet case, small cut-offs do not affect the momentum 
range in which first order perturbation theory works, al- 
though due to numerical limitations, the cut-off cannot 
be reliably reduced below r c = 0.3 fm. However, there 
are reasons for keeping the cut-off in the proposed win- 
dow, such as avoiding unphysical deeply bound states 
in the leading order amplitudes (the first one appears 
at r c = 0.45 fm), or an excessive D-state probability in 
the deuteron, yielding poor convergence in nuclear mat- 
ter calculations [62} . Larger cut-offs, of the order of f fm 
and above, are also disfavoured as they lead to a worse 
description of the e\ mixing angle for momenta above 
300 MeV, similar to the one obtained in the N 2 LO Wein- 
berg calculation of Ref. 17]. The proposed cut-off range 
avoids the previous problems and, due to the stronger 
long range distortion provided by the tensor component 
of OPE, generate error bands which decrease in size order 
by order. 

In this regard, it is interesting to notice the opposite 
cut-off preferences of the singlet and triplet channels. 
The mismatch in the preferred cut-off windows is a re- 
flection of the different physics at play in these waves. 
In the singlet, all pion exchanges are perturbative and 
the iteration of OPE is merely a short-cut to avoid the 
computation of higher order perturbations, while in the 
triplet tensor OPE really needs to be iterated. Different 
power countings require different cut-off windows. In this 
sense, large cut-off values worsen the convergence of the 
triplet: the OPE tensor force starts to behave perturba- 
tively, even if fully iterated in the Schrodinger equation. 
This entails a change in the counting of the triplet chan- 
nel from the modified Weinberg scheme of Nogga, Tim- 
mermans and van Kolck [l{| eventually to KSW [3^, [34| , 
thus reducing the convergence of the theory as the cut-off 
is increased. The exact point at which the change takes 
place is difficult to determine, but probably lies above 



8 This is to be compared with the singlet channel, where ~ 1 

and each subtraction adds an additional r 2 suppression. If the 
singlet LO potential had behaved as expected by power count- 
ing, i.e. 1/r 3 , it would have only needed two counterterms at 
NLO/N 2 LO, following Weinberg's counting. 
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FIG. 2. (Color online) Phase shifts (nuclear bar) for the 
3 Si - 3 £>i coupled channel, with the LO piece (OPE) fully 
iterated and the NLO and N 2 LO pieces (chiral TPE) treated 
perturbatively. The LO counterterm is fixed to reproduce the 
triplet scattering length ao,t = 5.419 fm. The error bands 
and the fitting range are the same as in Fig. fT}. The light 
blue band corresponds to the N 2 LO results of Ref. [TtJ in 
the standard Weinberg counting. The dashed dark blue line 
represents the N 2 LO results for r c — 0.3 fm. 



r c > 1/m-n = 1.4 fm. Of course, the fact that the ei mix- 
ing angle is the phase which starts to feel the problem 
earlier is not a surprise, as it depends on delicate cancel- 
lations between short and long range effects. The singlet 
channel, on the contrary, does not have any problem with 
larger values of the cut-off as the power counting is not 
changed: the use of larger cut-off values only entails a 
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rearrangement of the short range physics to account for 
those parts of the pion tail which have been ignored, but 
the assumption that all long range interactions are per- 
turbative remains unchanged. 



IV. THE ROLE OF THE CUT-OFF 

In the previous calculations we have taken a very prag- 
matic point of view with regard to the cut-off: we have 
chosen the cut-off range r c — 0.6 — 0.9 fm in order to im- 
prove the convergence of the theory and the description 
of the phase shifts up to N 2 LO. Of course, the choice of 
this range depends on a compromise between the specific 
requirements of the singlet and triplet channels, as ex- 
plained in the previous sections. The important point is 
however that the proposed cut-off window generates lead- 
ing order phase shifts which do not differ too much from 
the Nijmegen ones, an arrangement which minimizes the 
size of the subleading order corrections and, as a conse- 
quence, enhances the convergence of the theory. 

This criterion basically coincides with the interesting 
cut-off philosophy of Beane, Kaplan and Vuorinen [63j j . 
in which the cut-off is merely a parameter controlling the 
convergence rate of the theory. The underlying idea be- 
hind this interpretation is an analogy with the role of 
renormalization scale dependence in QCD (see for exam- 
ple [6J, |65|). A similar rationale can be provided by the 
observation that the full scattering amplitude, computed 
at all orders, is cut-off independent as a consequence of 
having an infinite number of counterterms. In this re- 
gard, cut-off dependence is just an artifact of finite order 
approximations which can be avoided by the careful se- 
lection of a cut-off window for which the particular power 
counting under consideration is realized. 

However, for this interpretation to be complete within 
an EFT context it is necessary to determine first some 
formal aspects of the theory, such as the expansion pa- 
rameter and the cut-off and momentum ranges for which 
the perturbative expansion converges. The knowledge of 
the expansion parameter is fundamental in order to be 
able to make rigorous error estimations of the results and 
to check the suitability of the selected cut-off window. 
On the other hand, the determination of the range of 
applicability of the EFT is necessary for avoiding power 
counting abuse, that is, claiming as legitimate the acci- 
dental description of data beyond the possibilities of the 
EFT under consideration. 

Due to the mostly numerical nature of the present 
investigation, it is not clear how to extract the expan- 
sion parameter. However, the deconstruction of the So 
singlet channel of Ref. [46] might provide some valu- 
able clues regarding this important aspect of the theory. 
The energy dependence of the short range physics in this 
channel suggests a breakdown scale of A 0iS ~ 270 MeV. 
This estimation translates into an expansion parameter 
of to w /Ao jS ~ 0.5, a value which is compatible with the 
conjectured equivalence of the present approach with the 



KSW counting in the singlet channel. For the 3 Si — 3 Di 
triplet channel there is no deconstruction yet which might 
provide a preliminary estimation of the breakdown scale, 
but if we assume the deconstruction of the p-wave un- 
coupled triplets [45| to hold in the 3 Si — 3 Di coupled 
triplet, we obtain Ao,t — 340 MeV 0. The related expan- 
sion parameter would be m ff /Ao.t — 0.4, a value which is 
compatible with the observation that the convergence in 
the triplet is better than in the singlet. 



V. RELATION TO OTHER APPROACHES 

In the present work we determine the power counting 
of the counterterms by requiring the renormalizability of 
the perturbative corrections to the scattering amplitude, 
where by renormalizability it is understood the elimina- 
tion of all negative powers of the coordinate space cut-off 
r c . There is still a residual cut-off dependence which 
is nominally of higher order, meaning that perturbative 
renormalizability implies the cut-off independence of the 
scattering amplitude at the order considered. 

This is very similar to the RG approach of Birse [58| , 
where the relative scaling (i.e. the power counting) of the 
counterterms is determined by requiring the cut-off inde- 
pendence of the scattering amplitude. Of course, exact 
cut-off independence is only achieved at infinite order. 
Finite order truncations will lead to a residual cut-off 
dependence involving positive powers of the cut-off r c , 
but the renormalizability of the amplitudes is guaran- 
teed. Therefore, it is not surprising a great degree of 
agreement between these two approaches. 

This expectation is realized in the singlet channel, 
where renormalization group analysis 



42| and decon- 



struction [46( | are equivalent to perturbative renormal- 
izability. For the triplet channel the situation is mixed: 
in the case of the t\ mixing angle and the 3 Z>i phase, the 
observation that two counterterms are needed to renor- 
malize each of these phases is compatible with the de- 
construction of the p- and d-wave uncoupled triplets of 
Ref. (4f|. However, the RGA of Ref. [42[ predicts one 
additional counterterm for the 3 Si phase which should 
appear at order Q 5 / 2 . This counterterm is not needed by 
perturbative renormalizability. 

The previous discrepancy is surprising: we are mak- 
ing the same assumptions as Ref. [42[ regarding which 
pieces of the interaction to iterate, yet the resulting power 
countings are slightly different. However, this is not new: 



9 This corresponds to a laboratory energy of 250 MeV, above 
which the short range interaction in the 3 Po, 3 Pi and 3 Z>2 waves 
cannot be reliably described by two counterterms l4, r | . On the 
other hand, the assumption that the p-waves yield a good ap- 
proximation for the breakdown scale of the s-waves is not unrea- 
sonable if we take into account that the deconstruction of the 
1 Pi wave [4^| basically suggests the same estimation as the 1 So 
wave |4g|. 
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the non-perturbative renormalizability of the OPE po- 
tential dictates that each attractive triplet requires one 
counterterm, while repulsive triplets do not. On the con- 
trary, the RGA of Ref. 42] makes no distinction for the 
power counting of attractive and repulsive triplets. The 
paradigmatic example is given by the 3 Po (attractive) 
and 3 Pi (repulsive) waves. As happened in the peripheral 
waves, the inconsistency can be circumvented in terms of 
the perturbative analysis of tensor OPE of Ref. [42[ : for 
the 3 Po wave the perturbative treatment of OPE is ex- 
pected to fail already at k ~ 200 MeV, while for the 3 Pi 
wave this limit is extended up to k ~ 400 MeV. There- 
fore, in the range of momenta of interest for nuclear EFT 
the 3 Pi wave can in principle be described in terms of the 
original Weinberg counting F^l . 

For the 3 Si phase the causes of the disagreement are 
to be found in the naive extrapolation of the idea of triv- 
ial and non-trivial fixed points to a problem where these 
concepts may not be applicable. The relevant observa- 
tion in this context is that attractive singular potentials 
do not have a unique solution [2(| HJ, HH : the value of 
the scattering length oscillates indefinitely as the cut-off 
varies, a situation which is solved by the inclusion of a 
counterterm, stabilizing the solution. In this regard, for 
an attractive singular interaction all values of the scat- 
tering length are equally fine-tuned, implying that the 
distinction between trivial and non-trivial fixed points is 
artificial in this case. As analyzed in Ref. the RG 
evolution of attractive singular potentials is driven in the 
infrared limit towards an oscillatory attractor-type so- 
lution resembling a limit cycle. However, the attractor- 
type solution does not have the discrete scaling properties 
of limit cycles (see [4(| for details) . 

Alternatively, the previous observations can also be un- 
derstood in terms of the behaviour of the squared reduced 
wave functions at short distances. For regular poten- 
tials there are two possible behaviours, the regular one, 
u(r c )| 2 ~ r 2 , which can be identified with "natural" sys- 
tems, and the irregular one, |w(r c )| 2 ~ 1, which describes 
systems with unnaturally large scattering lengths. On 
the contrary, for attractive singular potential, the wave 
function always behaves as |u(r c )| 2 ~ r 3 / 2 (times an os- 



10 One could think of extending this argument to the 3 D\ phase, 
which is usually well reproduced in perturbation theory [id . 
However, taking into account the coupled channel nature of the 
3 Z?l phase, it is probably inconsistent to treat tensor OPE per- 
turbatively in the d-wave channel but not in the s-wave channel. 

11 Notice however that Ref. [4fj| uses a different language than 
Birse's RGA [42j : what is called ultraviolet (long range) limit 
in 40] corresponds to the infrared (short range) limit of [42l |. If 
we call the light and heavy scales m; and , Ref. |4(I is taking 
m;r c — ¥ 0, while Ref. [42l ] assumes 3> l/r mi, or equiv- 
alently m^r e — > oo, m;r c — » and mj/mi, — > 0. Contrary to 
Ref. [42l |. Ref. [40l | does not analyze the power counting of the 
short range operators but rather concentrates on issues such as 
the cut-off dependence of observables and the fixed points, limit 
cycles and attractors which result from the RG flow of regular 
and singular potentials. 



cillatory factor), independently of the value of the scat- 
tering length. That is, there is no additional short range 
enhancement due to large scattering lengths. In this re- 
gard, we should not expect the existence of two different 
kinds of fixed points in the RG flow of attractive singular 
interactions. The previous observations indicate that for 
attractive singular potentials (i) the Co counterterm must 
be of order as required by non-perturbative renor- 

malizability and (ii) the first perturbation to the Co coun- 
terterms is of order Q -1 / 2 as expected from the behaviour 
of the squared wave function, meaning that the attractor 
is a stable solution of the RG flow. Consequently, the cor- 
rect RG analysis for channels with an attractive tensor 
force is the one termed "trivial" in Ref. [42[ , conveniently 
modified to incorporate the previous observation about 
the Co operator. 

A recent work which is also relevant for the present 
approach is the toy model proposed by Epclbaum and 
Gegelia to address the role of regularization and renor- 
malization in EFT 38]. In this work, the authors con- 
sider a two-body potential problem which shares many of 
the features of nuclear EFT, like the existence of a sep- 
aration of scales or the possibility of expanding the long 
range interaction in terms of a power counting. The con- 
clusions of the analysis of Epelbaum and Gegelia support 
most of the assumptions usually invoked in the Weinberg 
scheme, namely that naive dimensional analysis provides 
a good enough power counting and the ideal value of the 
cut-off should be chosen of the order of the hard scale 
of the problem. In addition, if the cut-off is taken much 
beyond the hard scale, the non-perturbative renormaliza- 
tion procedure may break the assumptions made in the 
first place by the power counting, a phenomenon which 
Epelbaum and Gegelia call "peratization" . 

The lessons derived from a specific toy model may be 
however of limited significance. In particular, there is 
an essential feature of the chiral expansion which is not 
reproduced in the previous model, namely the appear- 
ance of singular interactions at leading and subleading 
orders. Contrary to the expectations of Epelbaum and 
Gegelia, the presence of singular potentials implies that 
(i) non-perturbative power counting will break down at 
cut-offs much softer than expected and that (ii) devi- 
ations from naive dimensional analysis may eventually 
happen. These aspects have been probably overlooked in 
the previous analysis due to the very good properties of 
the toy model: subleading contributions to the toy po- 
tential are only mildly divergent and, in addition, they 
are always suppressed by the expected ratio of low energy 
versus high energy scales. On the contrary, the sublead- 
ing pieces of the chiral NN potential can receive unex- 
pectedly large contributions from light degrees of free- 
dom which have not been explicitly taken into account, 
like the A resonance. It is not surprising therefore that a 
toy model incorporating many of our naive expectations 
about EFT turns out to confirm them. 

However, as far as we limit ourselves to soft enough 
cut-offs, the conclusions of Epelbaum and Gegelia re- 
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garding naive dimensional analysis (i.e. Weinberg count- 
ing) are likely to hold. This observation is realized in 
the work of Shukla et al. [48| which, much in the spirit 
of deconstruction, analyzes the short distance physics of 
the singlet channel with the chiral NN potential up 
to N 2 LO. The authors observe that in the cut-off re- 
gion r c = 1.0 — 1.8 fm two counterterms are enough to 
parametrize the short range physics, a finding consis- 
tent with the idea that the Weinberg counting is bet- 
ter realized for soft values of the cut-off. A particularly 
interesting aspect of the previous work is the reanal- 
ysis of the short range physics for perturbative chiral 
TPE. For the cut-off range r c = 1.4 — 1.8 fm the ex- 
tracted short range physics can be accurately approxi- 
mated by first order perturbative TPE, while for the re- 
gion r c = 1.0 — 1.4 fm one needs to go to second and third 
order in the perturbative series in order to reproduce the 
non-perturbative results, although there is still conver- 
gence. In the softer cut-off range the Weinberg scheme is 
perfectly realized as a perturbative power counting. For 
the harder cut-off range, Weinberg is still a consistent 
(non-perturbative) power counting scheme, as sublead- 
ing order corrections are smaller than leading order ones. 
The efforts of Ref . [48j probably represent the best way to 
analyze the merits of the Weinberg counting in realistic 
cases. The extension to other partial waves, in particular 
the triplet, would be very welcomed. 

If the cut-off is decreased below R = 1.0 fm, the au- 
thors of Ref. [48[ observe that the contributions from sub- 
leading TPE start to grow uncontrollably, signalling the 
breakdown of the Weinberg counting. Below this cut-off, 
power counting is likely to be lost in non-perturbative 
calculations, as loop contributions from the subleading 
pieces will eventually dominate the amplitudes. The pre- 
vious breakdown scale is however uncomfortably soft: 
using the equivalence A = 7r/2r c [261 ] . i?o naively cor- 
responding to a (sharp) momentum cut-off of Ao ~ 
310 MeV. Most Weinberg calculations use momentum 
space cut-offs of the order of A ~ 0.5 GeV, which may be 
hard enough as to peratize the amplitudes. As suggested 
in Ref. [49J, this may be already happening in the 1 So 
singlet channel for A = 400 MeV at N 2 LO. These ob- 
servations do not imply however that Weinberg counting 
is not useful, only that it should be employed within its 
specific range of applicability. In this respect, the most 
interesting feature of perturbative treatments is that they 
are guaranteed to respect the power counting indepen- 
dently of the value of the cut-off, precluding from the 
start the possibility of any power counting inconsistency. 

A recent work which is also relevant for the discussion 
is the new KSW expansion of Beane, Kaplan and Vuori- 
nen [(33[ , which challenges one of the key premises of the 
present approach, namely that OPE should be fully it- 
erated in the triplet, by constructing a viable nuclear 
EFT in which all pion exchanges are treated as pertur- 
bations. In this work the convergence problems of the 
original KSW counting [42|, [55[ are alleviated by the ex- 
change of a fictitious meson of mass A which regulates 



the 1/r 3 singularity of the tensor force at short distances. 
For the optimum value of the regulator (A = 750 MeV), 
the expansion apparently converges up to order Q, al- 
beit slowly. At this order, the results of Ref. [63[ for the 
3 Si and 3 Di phases compare well with the LO results of 
the present approach. However, the order Q results for 
the ei mixing angle is clearly worse than our LO com- 
putation and it does not seem to converge for k > m^. 
This may be a good indicator that the tensor force re- 
ally needs to be iterated, as the t\ mixing angle is very 
sensitive to large cancellations between long and short 
range physics. In any case, a serious comparison of the 
present approach with the proposal of Beane, Kaplan and 
Vuorinen [63[ requires (i) the extension of the previous 
results beyond order Q and (ii) the consideration of the 
3 Po phase which according to Nogga, Timmermans and 
van Kolck (l9| also demands the non-perturbative inclu- 
sion of tensor OPE. 

The observation that OPE is perturbative in the sin- 
glet and non-perturbative in the triplet is closely related 
with the proposal of Beane, Bcdaquc, Savage and van 
Kolck (BBSvK) [54j, which suggested the iteration of 
those pieces of the (leading order) chiral NN potential 
which survive in the chiral limit (that is, tensor OPE). 
This prescription is theorized to generate a convergent 
expansion of the scattering amplitudes around the chi- 
ral limit, therefore providing a consistent EFT expansion 
for two-nucleon systems. The existence of a deeper rela- 
tionship with the present approach remains to be seen. 
However, the consideration of the subleading orders of 
the potential can break the correspondence, as there are 
pieces of these contributions to the potential which sur- 
vive in the chiral limit and which are strong enough as 
to be iterated, particularly in the singlet channel. In this 
regard, the BBSvK scheme might provide a justification 
for the iteration of chiral van der Waals forces. 



VI. CONCLUSIONS 

The present approach determines the momentum and 
cut-off range for which chiral TPE behaves perturba- 
tively when renormalizability is imposed. The use of 
small cut-offs is straightforward, but reduces the range 
of applicability of the theory in the singlet channel. The 
calculations turn out to confirm the viability of the count- 
ing proposal of Nogga, Timmermans and van Kolck [l9| , 
and corroborate to a large extent the related RGA by 
Birse [42] , which predicted the power counting of the 
short range operators. There are some minor discrepan- 
cies however between perturbative renormalizability and 
RGA in the triplet channel, specifically for the 3 Si phase, 
which are understood, suggesting minor modifications 
and possible improvements to the RGA of [12]. How- 
ever there are some formal aspects of the present EFT 
formulation which need to be elucidated, like the role of 
the cut-off, the determination of the expansion parame- 
ter or the range of applicability of perturbative TPE. In 
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this regard, the deconstruction approach of Refs. [45H47j 
is able to provide some interesting clues and preliminary 
answers. Of course, a complete evaluation of the renor- 
malized perturbative treatment of chiral TPE should also 
include the calculation of the p- and d-wave phase shifts 
and the deuteron properties. The present analysis paves 
the way for such computations, which we leave for future 
works. 
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Appendix A: Derivation of the DWBA for the 
Singlet Channel 

In this appendix, we derive the DBWA formulae used 
along the present paper. We start by considering a po- 
tential which can be decomposed as a zeroth order ap- 
proximation and a perturbation 

V(r) = V m (r) + V il \r), (Al) 

and the related reduced Schrodinger equations for the 



zeroth order and full reduced wave functions, u k 

Uk 



(0) 



and 



+ 2 M [V^ +V^]u k = k 2 



(A2) 

-U k - + 2(1 [V + V -''[ life =K'Uk, (A3) 

where (i is the reduced mass of the system. The full 
reduced wave function can be perturbatively expanded 
as 

(r) = uf{r) + 4 1} (r) + 0{{V^f) , (A4) 

where, for the purposes of this work, it is enough to con- 
sider first order perturbation theory only. 

In order to obtain the DWBA expressions we begin by 
(i) multiplying the zeroth order Schrodinger equation by 
the full solution u k , and (ii) the full Schrodinger equation 
by the zeroth order solution . Then we compute the 
difference between (i) and (ii), yielding 



,(°), 



,(o)' 



u k )' = 2(iV [1 \ r ) u k °\r)u k (r). 



The expression above can be integrated to obtain the 
Wronskian identity 



W(u k °\u k ) 



2(i 



drV^(r) 4 0) (r) u k {r) . 



(A6) 



where W(f, g) = f(r)g'(r) — f'{r)g(r) is the Wronskian, 
and r c and R are respectively the ultraviolet and infrared 
cutoffs. The infrared cutoff R can be eliminated by tak- 
ing into account the long distance behaviour of the u k °^ 
and u k reduced wave functions, which is given by 



u k {r) 



1 



sin (kr + 6(°1) 



AW(k) sin(5( ) 
1 sm(kr + S) 



A(k) 



sin S 



(A7) 
(A8) 



where A^ (fc) and A(k) are normalization factors which 
ensure an energy independent normalization of the re- 
duced wave functions at the cut-off radius. With the 
previous wave functions the Wronskian can be evaluated 
at R — > oo, resulting in 



W{ul 



(0) ,,(0+1) 



k sin(5-<5(°)) 
R ~ ~ A^A sin<5 sintfW 

Therefore, we arrive at the following expression 

k sin (6-6^) 
AW A sin 6 sin + -'^ ~ 

poo 

-2(i / drV^{r)uf{r)u k {r), 



(A9) 



(A10) 



where f(r c ) is just the Wronskian evaluated at r = r c , 
i.e. /(r c ) = W(u k \u k )\ ra , which does not depend on 
the momentum A: as a consequence of the energy inde- 
pendent normalization at r = r c . The perturbative ex- 
pansion of the previous formula can be obtained from the 
corresponding one of its components 

6(k) = 6 {0) (k) + S (1) (k) + 0((V {1) ) 2 ) , (All) 

u fc (r) = u k °\r) + u^(r) + 0((V^f) , (A12) 

A{k) = A (0) (k) + A w (k) + 0{{V {1) f) , (A13) 

f(r c ) = f {1) (rc) + 0((V^) 2 ), (A14) 

yielding the following DWBA formula for the phase shift 



A^ 2 sin5(0) 2 



-2(i 



drV {1) (r)u k °> (r), 



(A15) 



(A5) 



where the Wronskian term /W can be safely ignored in 
renormalized computations, as it vanishes once the first 
subtraction is done. 
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Appendix B: DWBA for the Triplet Channel 

In this appendix we present the perturbative distorted 
wave formulas for the phase shifts in the 3 <Si — 3 £>i triplet 
channel. For that, we express the phase shifts as the 
expansion 



^(fc;r c ) = 4 0) + 4 2) +4 3) 



,(3) 



-0(Q A ), 
- 0{Q A ) , 
0(Q 4 ), 



(Bl) 
(B2) 
(B3) 



e(k;r c ) 

where we have chosen the eigen parametrization of the 
phase shifts fo6| because in this parametrization the 
DWBA formulas take their simplest form. The expan- 
sion of the nuclear bar phase shifts [(57| can be obtained 
from the previous one by reexpanding the relationships 

Si+6 3 = S a +8p i (B4) 
sin (5 1 -5 2 ) = ^! ) (B5) 



sm(S a - Sp) 



tan 2e 
sin 2e 

sin2e ' 



according to the counting. The LO phase shifts, Sa \ 5^ 

and e(°\ are obtained by solving non-perturbatively the 
OPE potential with one counterterm, which is used for 
fixing the triplet scattering length to at = 5.419 fm. The 
exact procedure is explained in Ref . [6l[ . The expressions 
for the perturbative corrections to the LO phase shifts are 
the following 



sin 2 J< 0) 
SP{k;r c ) 



2fi 
~k 

2// 



A^\k)I^{k;r c ), 



~ 1.5 ^0^ (ty IppityTc) , 



sin 



(0) 





(B7) 



(B8) 



2 M Af\k)A?>(k) , v) 



k 3 cot4 0) -cot<5i 0) 



where the perturbative integrals I aa , 
defined as 



I$(k;r c )= dr 



Wire), 

(B9) 
and 7/3/3 are 



*#Hr)«S2(r)«jS (r) + 



(BIO) 



with p,<T — a, (3. As in the singlet case, \i represents the 



reduced mass of the system, u 



(0) 
k,a(0) 



and w 



(0) 

k, a {0) 



are the 



leading order s- and d-wave reduced wave functions for 
the a((3) scattering states in an energy independent nor- 
malization at the origin / cut-off radius, and Aa^ and 
Ag are the normalization factors which ensure the pre- 
vious condition. The asymptotic normalization of the a 



and P scattering states is taken to be 

4 0) «C(r) -> coseW (cottfW j (kr) - y (kr)) , 
4 0) w^Jr) -> sineW (cot j 2 (fcr) - fc(fcr)) , 



(Bll) 

fc2 -4 0) -»• - sine(0) ( cot 4 0) J'o( fcr ) - , 

fc 2 ^ 0) ^°J(r) -> cose (0) (cot4 0) i2(fcr) -fo(fcr)), 

(B12) 

where = xji(x), yi(x) = xy t (x), with j/(o;) and 

yz(x) the spherical Bessel functions. Due to the energy 
independent normalization of the wave functions at the 
cut-off radius, they can be expanded at short distances 

as [m 

OO 

4!L)W=E4 „U)« fc2 "' (Bi3) 



(B6) 

where the behaviour is given by 



i (0) 

l 2n,a(0) 



(r) 



r 3/4+5n/2 



(B14) 



with f(x) some combination of sinx, cosx and e~^ x . 
The length scale a is related to the strength of the tensor 
force. In principle the general solution of the Schrodinger 
equation for the tensor OPE potential also admits an 
e +V2x com p nent which would destroy the renormaliz- 
ability of the theory, as it generates divergences which 
cannot be absorbed by any finite number of countcrt- 
erms. The previous component does not appear however 
if the LO wave functions have been properly renormal- 
ized. Therefore, what is essential is the power law be- 
haviour of the wave functions, which dictates the diver- 
gence structure of the perturbative integrals 1$) 



/W(fc;r c )=7^ ff (r c ) + fc 2 /^(r c ) 



0,per 
+ lR,pa(k'> r c) 



(B15) 



with io^po- and jj" ' the divergent pieces of the integral 
and Ift pcr the regular piece. We can regularize the inte- 
gral I^pa by including two free parameters 

/&> (k; r c ) = Ag CT + A^ CT fc 2 + I&> (k; r c ) , (B16) 

which are to be fitted to the scattering data of the cor- 
responding phase. The previous procedure yields a total 
of six counterterms for regularizing the NLO and N 2 LO 
phase shifts. As in the singlet case, the finite piece of 
one of these parameters (Aq is redundant as it only 
affects the triplet scattering length at , which was already 
fixed at leading order. In other words, six pieces of data 
are enough to determine the NLO/N 2 LO phase shifts in 
the triplet. 
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